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We present the first measurement of the mass of the top quark in a sample of tt (Mbbqq events 
(where t = e, ^) selected by identifying jets containing a niuon candidate from the seniileptonic 
decay of heavy-flavor hadrons (soft muon 6-tagging). The pp collision data used corresponds to an 
integrated luminosity of 2 fb^^ and was collected by the CDF II detector at the Fermilab Tevatron. 
The measurement is based on a novel technique exploiting the invariant mass of a subset of the 
decay particles, specifically the lepton from the W boson of the t Wb decay, and the muon from 
a semileptonic b decay. We fit template histograms, derived from simulation of tt events and a 
modeling of the background, to the mass distribution observed in the data and measure a top quark 
mass of 180.5 ± 12.0(stat.) ± 3.6(syst.) GeV/c^, consistent with the current world average. 
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PACS numbers: 14.65.Ha, 12.15.Ff 



A massive top quark plays an important role in the 
standard model (SM). The mass of the top quark {mt) 
enters electroweak (EW) precision observables as an in- 
put parameter via quantum effects, i.e. loop corrections, 
and its large numerical value gives rise to sizable correc- 
tions that behave as powers of mj f]\. For example, in the 
theoretical prediction of the W boson mass (mw) within 
the SM, when these corrections are combined with the 
logarithmic dependence on the mass of the postulated 
Higgs boson (m/f), a relationship emerges that provides 
a constraint on mn from experimental determinations of 
mw and rrit [1]. Indeed, the strong dependence of the 
SM radiative corrections on rrit made it possible to pre- 
dict the value of mt Q prior to its experimental determi- 
nation [1, @. Thus, a precision value of Wt is crucial for 
constraining SM parameters, for high-sensitivity searches 
for effects of new physics and for stringent consistency 
tests of models beyond the SM (e.g. supersymmetry). 
Furthermore, independent measurements of mt in all fi- 
nal states of tt decay provide an important consistency 
check of the top quark sector of the SM, and might reveal 
new physics with top-like signatures. 

Significant progress has been made recently in reduc- 
ing the uncertainty in measurements of mt and in devis- 
ing alternative and independent techniques. The current 
best single measurement is determined by reconstructing 
the full decay chain and computing the invariant mass 
of the decay products in tt iPbbqq events, and yields 
mt = 172.1 ± 1.6 GeV/c2 d, Q. However, this and aU 
the most precise of the current techniques are limited by 
the common systematic uncertainty in the calorimeter 
jet energy calibration (jet energy scale, JES). To pro- 
vide independent measurements, several techniques with 
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minimal dependence on the JES have been proposed. For 
example, the flight distance of the 6-hadron from the top 
decay can be used to infer the mass of the top quark 
Q, but this method also requires precision track recon- 
struction to determine the decay length. A proposal has 
been made for exploiting the correlation between mt 
and the invariant mass of the system composed of a J/^p 
(from the decay of a 5 hadron) and the Icpton from the 
W decay. The advantage is a stronger correlation of this 
system-mass with mt than that of individual decay prod- 
ucts of the top quark, and thus a better sensitivity to the 
top quark mass, but the overall branching ratio for this 
final state is only O(10~^). 

We present the first measurement of the mass of the top 
quark in a sample of ti (Dbbqq events (where i = e, fj,) 
selected by identifying 6-jets with a candidate muon from 
semileptonic decay of heavy-fiavor hadrons. We have de- 
veloped a novel technique that exploits the invariant mass 
of the lepton from the W boson of the t Wb decay, 
and the muon from a semileptonic b decay. The selection 
method is complementary to that taking advantage of the 
long lifetime of 6-hadrons through the presence of a decay 
vertex displaced from the primary interaction. Since only 
^ 50% of the sample of tt candidates with a semileptonic 
b decay overlaps the top samples selected by the identifi- 
cation of a displaced vertex, and a still smaller fraction is 
in common with traditional samples that require all four 
jets for the mass reconstruction, our technique provides 
an essentially independent measurement of mt from these 
data. Moreover, our observable is largely independent 
of the JES, because the calorimeter information is used 
solely for the selection of event candidates, and therefore 
the result can add a significant amount of information 
when averaged with those from other measurements. In- 
cluding sequential decays of charm, the branching frac- 
tion for b — ^ fii^X ~ 20% is sizable and since this 
technique does not require precision secondary vertex re- 
construction to suppress backgrounds, it could be an at- 
tractive option for the early phase of experiments at the 
Large Hadron Collider (LHC). Finally, the observable has 
a higher correlation to the top quark mass than the mo- 
mentum of the lepton from the W decay alone. A partial 
reduction in sensitivity will arise from b-W mis-pairing, 
when the lepton from the W decay and the muon from 
the b semileptonic decay do not originate from the same 
top quark. 

Top quarks are produced at the Tevatron proton- 
antiproton collider predominantly in pairs of t and t, and 
are identified by the SM decay t Wb, providing a fi- 
nal state that includes two W bosons and two bottom 
quarks. W^s are identified through their decay to lep- 
tons or quarks. Quarks hadronize and are observed as 
jets of charged and neutral particles. The CDF II de- 
tector is described in detail elsewhere [l^. The compo- 
nents relevant to this analysis include the central outer 
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tracker (COT), the central electromagnetic and hadronic 
calorimeters, the central muon detectors and the lumi- 
nosity counters. The data sample, produced in pp col- 
lisions at -^i = 1.96 TeV during Run II of the Fermi- 
lab Tevatron, was collected between March 2002 and 
May 2007 and corresponds to an integrated luminosity 
of 2.0 ±0.1 fb^^. We select events where one of the 
W bosons decays to an isolated electron (muon) carry- 
ing large transverse energy {Et) (momentum (pt)) |11| 
with respect to the beam line, plus a neutrino. We refer 
to these high-p-r electrons or muons as primary leptons 
(PL). The neutrino escapes the detector causing an im- 
balance of total transverse energy vector, referred to as 
missing Et [tr)- The other W boson in the event de- 
cays hadronically to a pair of quarks. We take advan- 
tage of the semileptonic decay of B hadrons by search- 
ing for muons within final-state jets (soft-lepton tagging, 
or SLT), in order to identify those jets that result from 
hadronization of the bottom quarks. 

The event selection starts with an inclusive lepton trig- 
ger requiring an electron (muon) with Et > 18 GeV 
(Pt > 18 GeV/c). Further selection requires that can- 
didate electron (muon) PLs are isolated and have Et > 
20 GeV (pt > 20 GeV/c) and It; | < 1.1. We define an 
isolation parameter, /, as the calorimeter transverse en- 
ergy in a cone of opening AR = \J {AriY + {Atpy — 0.4 
around the lepton (not including the lepton energy itself) 
divided by the electron Et or muon pt- We select iso- 
lated electrons (muons) by requiring / < 0.1. The event 
must have (t > 30 GeV, consistent with the presence 
of a neutrino from the W boson decay. Jets are iden- 
tified using a fixed-cone algorithm with a cone opening 
of AR — 0.4 and are constrained to originate from the 
pp collision vertex. Muons inside jets are identified by 
matching the tracks of the jet, as measured in the COT, 
with track segments in the muon detectors. Such a muon 
with Pt > 3 GeV/c and within AR < 0.6 of a jet axis 
is called an SLT/x [T3|. The probability of misidentify- 
ing a hadron as an SLT/i, denoted as the SLT/x mistag 
probability, is measured using a data sample of pions, 
kaons and protons from D* and A° decays. A Monte 
Carlo (MC) simulation of -flight flavor events is used 
to model the tt, K and p admixture in light-quark jets. 
The SLT/x mistag probability is parametrized as a func- 
tion of the track pT and ?/, and is seen to describe within 
±5% the number of false SLT^ tags in light flavor jets of 
QCD multijet and 7 -I- jet events. 

To reduce background from dimuon resonances and 
double-semileptonic B hadron decays, we remove events 
in which the PL muon and SLT/i are oppositely charged 
and have an invariant mass consistent with a Z, T 
or, irrespectively of the PL flavor, less than 5 GeV/c^. 
We further reject events as candidate radiative Drell- 
Yan and Z bosons if the tagged jet has an electromag- 
netic energy fraction above 0.8 and only one track with 
Pt > 1.0 GeV/c within a cone of AR = 0.4 about the jet 
axis The jet energies are corrected to account for varia- 
tions of the detector response in 77 and time, calorime- 



ter gain drifts, non linearity of calorimeter energy re- 
sponse, multiple pp interactions in an event and for en- 
ergy loss in un- instrumented regions [2l[. Finally, the 
sample is partitioned according to the number of jets with 
Et > 20 GeV and \ri\ < 2.0 in the event, and at least one 
jet is required to contain an SLT/i (defining the SLTfi- 
tagged W+n jets sample). The subset of W plus at least 
3 jets is the tt candidate sample, and to reduce back- 
ground from QCD production of W with multiple jets, 
we additionally require the total transverse scalar energy 
in the event {Ht [231) to be greater than 200 GeV. 

Standard model processes that result in the same sig- 
nature as the tt signal are backgrounds to this mea- 
surement. There are three dominant backgrounds: the 
largest one is mistags of -flight flavor events, and a 
smaller contribution is due to W boson in association 
with heavy flavor jets {Wbb, Wcc, Wc). Events without 
W bosons that pass the event selection are typically QCD 
multijet events where one jet has been reconstructed as a 
high-pT lepton, mismeasured jet energies produce appar- 
ent f^T ^^^d an additional jet contains an SLT/i. A frac- 
tion of these events is from bb and cc, where the candidate 
PL may result from a semileptonic decay of one of the 
fragmenting heavy quark and the SLT// from a semilep- 
tonic decay of the other. Other minor backgrounds that 
can mimic a W boson and an SLT/i signature include 
diboson {WW, ZZ, WZ), Drell-Yan^ tt, single top 
quark, and residual Drell-Yan— > /x/z events not removed 
by the dimuon resonance removal. The composition of 
the data sarnple used in this analysis has been studied ex- 
tensively in [I4I , where we have measured the production 
cross section for pp — > ttX, and is summarized in Table [D 
The W^-l-jets, QCD multijet and Drell-Yan background 
are determined using the data, while the remaining back- 
grounds are estimated from MC simulations. The W+1,2 
jets samples contain little tt events and have a com- 
position similar to the background of the tt candidate 
sample. The simulation of tt events is performed us- 
ing PYTHIA [l3| and HERWIG [ll|- The generators are 
used with the CTEQ5L [l^ parton distribution func- 
tions (PDF). Modehng of b and c hadron decay is pro- 
vided by EVTGEN [1^. Modeling of W-fjets production is 
performed using ALPGEN [i3|, coupled with pythia for 
the shower evolution and evtgen for the heavy flavor 
hadron decays. Diboson production {WW, ZZ, WZ) 
and Drell-Yan^ tt are determined using pythia. Drell- 
Yan^ /i/i -f jets events are modeled using ALPGEN while 
single top production is modeled with MADEVENT p^ . 
both with PYTHIA showering. The CDF II detector sim- 
ulation models the response of the detector to particles 
produced in pp collisions. The detector geometry used in 
the simulation is the same as that used for reconstruction 
of the collision data. Details of the CDF II simulation, 
based on the GEANtS package, can be found in (19l] . 

We compute the invariant mass {Mi^) between the PL 
and the SLT/i in the ti candidates sample. In rare cases 
where there is more than one SLT/i tag in the same jet, 
or more than one SLT/i tagged jet in the same event, we 
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TABLE L Composition of the SLT/i-tagged W + n jets candi- 
date sample [l^l • The Ht > 200 GeV requirement is released 
for events with fewer than 3 jets. 



W+1 ,2 jets sample 



Source 


W+1 jet 


W+2 jet 


W+ > 3 jets 


VP'+light flavor 


622±31 


226±12 


52.3±2.6 


ly+heavy flavor 


145±55 


66.6±25.2 


14.3±5.4 


QCD multijet 


91.9±16.5 


44.9±10.4 


6.9±1.5 


WW + WZ + ZZ 


3.8±0.4 


7.0±0.7 


1.9±0.3 


Drell-Yan^ tt 


2.6±0.6 


1.5±0.4 


0.6±0.3 


Drell-Yan— > fifj, 


6.0±1.2 


4.1±0.9 


0.8±0.5 


Single top 


4.4±0.4 


9.0±0.7 


2.7±0.2 


Total background 


876±54 


359±24 


79.5±5.3 


tt {att = 9.1 pb) 


3.5±0.2 


31.8±1.0 


168.5±5.3 


Data 


892 


384 


248 



use the SLT// candidate that has the best match between 
the COT track and the track segment in the muon detec- 
tors. No attempt is made to choose the correct pairing 
from the decay chain of the two top-quarks. The elec- 
tric charge of the SLT/x for instance is not an effective 
flavor selector due to abundant sequential b ^ c ^ fi 
decays. When the wrong pairing is chosen, there is still 
sensitivity to the top quark mass due to the boost of the 
SLT^ and the PL. The distribution of AIi^ is given by the 
contribution of ti and background events. For the back- 
ground, the Mi^ distribution of QCD multijet events is 
derived from the data themselves in the kinematic-region 
of / > 0.15, > 30 GeV, topologically close to the sig- 
nal region, while for other background sources we use MC 
simulation. We check the background model in W+1,2 
jet SLT^-tagged data events, a sample with a similar 
composition as the background to tt candidates. We find 
the predicted and observed distributions of M^^, (Figure 
[T|) to be in agreement with a p- value of 55%, as given by 
the Kolmogorov-Smirnov test. 

We construct a set of template histograms of the M^^ 
distribution using the background model and a simula- 
tion of events. The tt samples are generated with differ- 
ent top quark mass values in the range 150-195 GeV/c^, 
incrementing by steps of up to 0.5 GeV/c^, and the full 
Mi^ spectra are determined by adding the signal and ex- 
pected background histograms in the ratio shown in Ta- 
ble HI Figure [2] shows the mean value of the M^^ distribu- 
tions versus the input top quark mass, indicating a linear 
relationship between the two quantities. Also shown is 
< Mifj^ >:= 35.6 ± l.l(stat.) GeV/c^, measured in the 
data. We perform a binned-likelihood fit to the M^^ his- 
togram of the data, in 20 bins between 4-100 GeV/c^, 
with the binning and range chosen a priori appropriately 
to the size of the data sample. The likelihood is defined 
as: 



lnL(m() 



Nbins 

E 



^data 



In 



,TP 



{mt) 



tTP 



(1) 



• Data, 2 fb"' 

I I W+jets 

I I QCD multijet 

^ other 




10 20 30 40 50 60 70 80 90 100 

/W,„ [GeV/c^] 



FIG. 1: The predicted and observed M^^ distributions in the 
sample of W-f 1,2 jet SLT/^-tagged events. The predicted dis- 
tributions are stacked. 



each i-bin of the data and template histograms respec- 
tively, the total number of entries is nj^ = n'^ot'^, and 
nj^{mt)/ri^^ = Vi{mt) is the probability of the z-th bin, 
normalized such that J^i'^i — 1- The background nor- 
malization is fixed and its value is varied in the evaluation 
of the systematic uncertainty. A parabolic function is fit 
to the values of InL(mt) derived from each mass template, 
and the measured top quark mass is determined from the 
minimum of the likelihood function, while the statistical 
uncertainty is given by the range corresponding to an in- 
crease in the — InL of 0.5 units above the minimum. For 
each mass point within the full mass range, we gener- 
ate 5000 pseudoexperiments with the same sample size 
as that of the data and verify that the fitting procedure 
is unbiased and that the statistical uncertainty returned 
by the fits represents the 68% confidence level. From 248 
tt candidate events, we measure: 



mt 



180.5 ± 12.0(stat.) ± 3.6(syst.) GeV/c^ 



(2) 



where nf^^^ and nj^ (mt) are the number of entries in 



Figure[3]shows the AIi^ distribution of the data, the back- 
ground, and the templates corresponding to the best fit 
and the statistical uncertainty. 

The sources of systematic uncertainty that affect the 
measured value of the top quark mass are summarized 
in Table [Til The limited size of the tt samples simulated 
with different values of m^, input to the fitting procedure, 
yields an uncertainty of ±0.3 GeV/c^. Several compo- 
nents enter the uncertainty on the modeling of the back- 
ground. The uncertainty on the W+ heavy and light fia- 
vor normalizations yields an uncertainty of ±0.5 GeV/c^. 
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MC simulation 
Fit to simulation 
Data ±1 o, 2 fb^ 



150 160 170 180 190 

Top quark mass [GeV/(f] 



FIG. 2: The correlation between the mean value of the M^^ 
histograms from simulated tt and background samples, and 
the input nit. The continuous line shows a linear fit to the 
points. 



• Data, 2fb ' 

□ Fit to data 

(m= 180.5 GeV/c'') 

Fit +1 a (Stat.; 

■ ■ ■ Fit -1 a (Stat.) 

□ Background 




10 20 30 40 50 60 70 80 90 100 



Mi^[GeV/c^] 



FIG. 3: The distribution of invariant mass M^^ of the lep- 
ton from the W decay and the SLT/x, from a sample of 248 
candidate tt events with 79.5 background. 



The uncertainty on the shape of the PF+jets histogram 
is evaluated by varying the distribution, to within the 
statistical accuracy associated with the comparison in 
the W+1,2 jets sample between the data and the back- 
ground model, and yields an uncertainty of ±1.4 GeV/c^. 
The normalization of the QCD multijet background con- 
tributes ±0.8 GeV/c^ The shape of the QCD multijet 



distribution accounts for ±0.6 GeV/c^, as determined by 
replacing the nominal sample with dijet enriched data 
selected by / < 0.1 and < 15 GeV, and by varying 
the distribution according to its statistical uncertainty. 
The shift on the measured top quark mass due to the 
uncertainties on the remaining backgrounds is negligi- 
ble. The total uncertainty from background modeling is 
±1.9 GeV/c2. 

Monte Carlo modeling of the signal Algf^ distributions 
includes effects of PDFs, initial-state radiation (ISR), 
final-state radiation (FSR), and JES. The uncertainty 
due to the MC modeling of ti production and decay, in- 
cluding b fragmentation, is determined by comparing the 
simulation using pythia with that using herwig and 
gives Amt — ±2.1 GeV/c^. The PDF uncertainty is 
evaluated by adding in quadrature the contribution of 
four effects: variations of the PDFs according to the 20 
CTEQ eigenvectors [1^ , the difference between the stan- 
dard tt simulation using the CTEQ5L PDF and one de- 
rived using MRST98 [2^ in the default configuration or 
with two alternative choices for Og, and the variation of 
the contribution of gluon fusion in tt production between 
5 and 20%. The overall estimated uncertainty from PDF 
is ±1.0 GeV/c2._We vary both ISR and FSR simulta- 
neously in the tt Monte Carlo simulation, within con- 
straints set by studies of radiation in Drell-Yan events in 
the data, and assign a systematic uncertainty on mt of 
±1.3 GeV/c2. 

The jet reconstruction is used in this analysis only for 
the selection of event candidates and therefore the un- 
certainty on the calibration of the jet energies enters the 
measurement solely through the event selection, via the 
jet counting and the requirement. The uncertainty 
due to the JES is measured by shifting the energies of 
the jets in tt MC simulation by ±lcr of the JES [2ll 
and results in Arrit = ±0.3 GeV/c^. The uncertainty 
of ±1% on the difference between data and simulation 
of the PL energy and momentum scales gives an uncer- 
tainty of ±0.9 GeV/c^. The differences in the data ver- 
sus simulation for the SLT/i px spectrum depends on the 
6-quark fragmentation modeling and the momentum cal- 
ibration. In addition to the different fragmentation mod- 
els in HERWIG versus pythia, we consider comparisons of 
the data with MC simulation of the muon px spectra in 
B ^-D°X [Hi and bb niiX [H] which indicate an 
uncertainty on the muon px of ±0.8%, corresponding 
to Arrit = ±0.9 GeV/c^. The uncertainty on the px de- 
pendence of the SLT/i tagging efficiency yields a shift on 
the top quark mass of ±0.2 GeV/c^. Finally, a source of 
systematic uncertainty is due to the modeling of pile-up 
events from multiple pp interactions and it is estimated 
to affect the measured mass by < ±0.5 GeV/c^. 

In summary, we have performed the first measurement 
of the top quark mass in a sample of tt — > £Dbbqq events 
selected by identifying 6-jets with a muon candidate from 
the semileptonic decay of heavy- flavor hadrons. The re- 
sult, mt = 180.5 ± 12.0(stat.) ± 3.6(syst.) GeV/c^, is 
in agreement with the current world average value of 
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TABLE IL Summary of systematic uncertainties. 



Source 


Amt [GcV/c ] 


MC tt samples statistics 


±0.3 


Background 


±1.9 


tt production and decay model 


±2.1 


Parton distribution functions 


±1.0 


Initial- and final-state radiation 


±1.3 


Jet energy scale 


±0.3 


PL energy /momentum scale 


±0.9 


SLT^ momentum 


±0.9 


Pileup 


±0.5 


Total 


±3.6 



173.1±1.3 GeV/cP 0, providing a consistency check of 
the top quark sector with soft muon 6-tagged events. Our 
measurement technique exploits the correlation between 
the parent top quark mass and the invariant mass of the 
system composed of the lepton from the W decay and the 
muon from the semileptonic B decay. The uncertainty at 
present is dominated by the statistical component. The 
method has a minimal dependence on the jet energy cal- 
ibration, making it suitable for averaging the result with 
those from other techniques, and its dominant system- 



atic uncertainties are likely reducible, e.g. by improving 
the calibration of the leptons' pr to better than 1% with 
J/f/;, T and Z resonances, by using improved tuning for 
the MC modeling of tt production and decay, and with 
high statistics data samples for the background model. 
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